Methods and Apparatuses for Signaling with Geometric Constellations in a Raleigh Fading Channel

ABSTRACT

Communication systems are described that use signal constellations, which have unequally spaced (i.e. ‘geometrically’ shaped) points. In many embodiments, the communication systems use specific geometric constellations that are capacity optimized at a specific SNR, over the Raleigh fading channel. In addition, ranges within which the constellation points of a capacity optimized constellation can be perturbed and are still likely to achieve a given percentage of the optimal capacity increase compared to a constellation that maximizes d min , are also described. Capacity measures that are used in the selection of the location of constellation points include, but are not limited to, parallel decode (PD) capacity and joint capacity.

CROSS-REFERENCE TO RELATED APPLICATIONS

The present invention is a continuation of U.S. patent application Ser. No. 14/943,003 entitled “Methods and Apparatuses for Signaling With Geometric Constellations in a Raleigh Fading Channel” to Barsoum et al., filed Nov. 16, 2015, which application is a continuation of U.S. patent application Ser. No. 13/179,383 entitled “Methods and Apparatuses for Signaling With Geometric Constellations in a Raleigh Fading Channel” to Barsoum et al., filed Jul. 8, 2011 which issued as U.S. Pat No. 9,191,148, which application claims priority as a Continuation-In-Part to U.S. patent application Ser. No. 12/156,989 entitled “Design Methodology and Method and Apparatus for Signaling with Capacity Optimized Constellation”, which issued as U.S. Pat. No. 7,978,777 and claims priority to U.S. Provisional Application Ser. No. 60/933,319 entitled “New Constellations for Communications Signaling: Design Methodology and Method and Apparatus for the New Signaling Scheme” to Barsoum et al., filed Jun. 5, 2007. U.S. patent application Ser. No. 13/179,383 entitled “Methods and Apparatuses for Signaling With Geometric Constellations in a Raleigh Fading Channel” to Barsoum et al., filed Jul. 8, 2011 which issued as U.S. Pat. No. 9,191,148 also claims priority to U.S. Provisional Application Ser. No. 61/362,649 filed Jul. 8, 2010 entitled “Methods and Apparatuses for Signaling with Geometric Constellations in a Raleigh Fading Channel” to Barsoum et al. The disclosures of U.S. patent application Ser. Nos. 14/943,003, 13/179,383, 12/156,989 and U.S. Provisional Application Nos. 60/933,319 and 61/362,649 are expressly incorporated by reference herein in their entirety.

STATEMENT OF FEDERALLY SPONSORED RESEARCH

This invention was made with Government support under contract NAS7-03001 awarded by NASA. The Government has certain rights in this invention.

BACKGROUND

The present invention generally relates to bandwidth and/or power efficient digital transmission systems and more specifically to the use of unequally spaced constellations having increased capacity.

The term “constellation” is used to describe the possible symbols that can be transmitted by a typical digital communication system. A receiver attempts to detect the symbols that were transmitted by mapping a received signal to the constellation. The minimum distance (d_(min)) between constellation points is indicative of the capacity of a constellation at high signal-to-noise ratios (SNRs). Therefore, constellations used in many communication systems are designed to maximize d_(min). Increasing the dimensionality of a constellation allows larger minimum distance for constant constellation energy per dimension. Therefore, a number of multi-dimensional constellations with good minimum distance properties have been designed.

Communication systems have a theoretical maximum capacity, which is known as the Shannon limit. Many communication systems attempt to use codes to be able to transmit at a rate that is closer to the capacity of a communication channel. Significant coding gains have been achieved using coding techniques such as turbo codes and LDPC codes. The coding gains achievable using any coding technique are limited by the constellation of the communication system. The Shannon limit can be thought of as being based upon a theoretical constellation known as a Gaussian distribution, which is an infinite constellation where symbols at the center of the constellation are transmitted more frequently than symbols at the edge of the constellation. Practical constellations are finite and transmit symbols with equal likelihoods, and therefore have capacities that are less than the Gaussian capacity. The capacity of a constellation is thought to represent a limit on the gains that can be achieved using coding when using that constellation.

Prior attempts have been made to develop unequally spaced constellations. For example, a system has been proposed that uses unequally spaced constellations that are optimized to minimize the error rate of an uncoded system. Another proposed system uses a constellation with equiprobable but unequally spaced symbols in an attempt to mimic a Gaussian distribution.

Other approaches increase the dimensionality of a constellation or select a new symbol to be transmitted taking into consideration previously transmitted symbols. However, these constellation were still designed based on a minimum distance criteria.

SUMMARY OF THE INVENTION

Systems and methods are described for constructing a modulation such that the constrained capacity between a transmitter and a receiver approaches the Gaussian channel capacity limit first described by Shannon [ref Shannon 1948]. Traditional communications systems employ modulations that leave a significant gap to Shannon Gaussian capacity. The modulations of the present invention reduce, and in some cases, nearly eliminate this gap. The invention does not require specially designed coding mechanisms that tend to transmit some points of a modulation more frequently than others but rather provides a method for locating points (in a one or multiple dimensional space) in order to maximize capacity between the input and output of a bit or symbol mapper and demapper respectively. Practical application of the method allows systems to transmit data at a given rate for less power or to transmit data at a higher rate for the same amount of power.

One embodiment of the invention includes a transmitter configured to transmit signals to a receiver via a communication channel, where the transmitter, includes a coder configured to receive user bits and output encoded bits at an expanded output encoded bit rate, a mapper configured to map encoded bits to symbols in a symbol constellation, a modulator configured to generate a signal for transmission via the communication channel using symbols generated by the mapper, where the receiver, includes a demodulator configured to demodulate the received signal via the communication channel, a demapper configured to estimate likelihoods from the demodulated signal, and a decoder that is configured to estimate decoded bits from the likelihoods generated by the demapper. In addition, the symbol constellation is a PAM-8 symbol constellation having constellation points within at least one of the ranges specified in FIGS. 25-47.

In a further embodiment, the code is a Turbo code. In another embodiment, the code is a LDPC code.

In a still further embodiment, the constellation provides an increase in capacity over the Raleigh fading channel at a predetermined SNR that is at least 5% of the gain in capacity achieved by a constellation optimized for joint capacity at the predetermined SNR.

In still another embodiment, the constellation provides an increase in capacity over the Raleigh fading channel at a predetermined SNR that is at least 20% of the gain in capacity achieved by a constellation optimized for joint capacity at the predetermined SNR.

In a yet further embodiment, the constellation provides an increase in capacity at a predetermined SNR over the Raleigh fading channel that is at least 50% of the gain in capacity achieved by a constellation optimized for joint capacity at the predetermined SNR.

In yet another embodiment, the constellation provides an increase in capacity over the Raleigh fading channel at a predetermined SNR that is at least 90% of the gain in capacity achieved by a constellation optimized for joint capacity at the predetermined SNR.

In another embodiment again, the constellation provides an increase in capacity over the Raleigh fading channel at a predetermined SNR that is at least 100% of the gain in capacity achieved by a constellation optimized for joint capacity at the predetermined SNR.

In a further additional embodiment, the constellation provides an increase in capacity over the Raleigh fading channel at a predetermined SNR that is at least 5% of the gain in capacity achieved by a constellation optimized for PD capacity at the predetermined SNR.

In another additional embodiment, the constellation provides an increase in capacity over the Raleigh fading channel at a predetermined SNR that is at least 20% of the gain in capacity achieved by a constellation optimized for PD capacity at the predetermined SNR.

In a still yet further embodiment, the constellation provides an increase in capacity over the Raleigh fading channel at a predetermined SNR that is at least 50% of the gain in capacity achieved by a constellation optimized for PD capacity at the predetermined SNR.

In still yet another embodiment, the constellation provides an increase in capacity over the Raleigh fading channel at a predetermined SNR that is at least 90% of the gain in capacity achieved by a constellation optimized for PD capacity at the predetermined SNR.

In still another embodiment again, the constellation provides an increase in capacity over the Raleigh fading channel at a predetermined SNR that is at least 100% of the gain in capacity achieved by a constellation optimized for PD capacity at the predetermined SNR.

A still further additional embodiment includes a transmitter configured to transmit signals to a receiver via a communication channel, where the transmitter, includes a coder configured to receive user bits and output encoded bits at an expanded output encoded bit rate, a mapper configured to map encoded bits to symbols in a symbol constellation, a modulator configured to generate a signal for transmission via the communication channel using symbols generated by the mapper, where the receiver, includes a demodulator configured to demodulate the received signal via the communication channel, a demapper configured to estimate likelihoods from the demodulated signal, and a decoder that is configured to estimate decoded bits from the likelihoods generated by the demapper. In addition, the symbol constellation is a PAM-16 symbol constellation having constellation points within at least one of the ranges specified in FIGS. 48-71.

Still another additional embodiment includes a transmitter configured to transmit signals to a receiver via a communication channel, where the transmitter, includes a coder configured to receive user bits and output encoded bits at an expanded output encoded bit rate, a mapper configured to map encoded bits to symbols in a symbol constellation, a modulator configured to generate a signal for transmission via the communication channel using symbols generated by the mapper, where the receiver, includes a demodulator configured to demodulate the received signal via the communication channel, a demapper configured to estimate likelihoods from the demodulated signal, and a decoder that is configured to estimate decoded bits from the likelihoods generated by the demapper. In addition, the symbol constellation is a PAM-32 symbol constellation having constellation points within at least one of the ranges specified in FIGS. 72-95.

Another further embodiment includes a transmitter configured to transmit signals to a receiver via a communication channel, where the transmitter, includes a coder configured to receive user bits and output encoded bits at an expanded output encoded bit rate, a mapper configured to map encoded bits to symbols in a symbol constellation, a modulator configured to generate a signal for transmission via the communication channel using symbols generated by the mapper, where the receiver, includes a demodulator configured to demodulate the received signal via the communication channel, a demapper configured to estimate likelihoods from the demodulated signal, and a decoder that is configured to estimate decoded bits from the likelihoods generated by the demapper. In addition, the symbol constellation is a N-Dimensional symbol constellation, where the constellation points in at least one dimension are within at least one of the ranges specified in FIGS. 25-95.

BRIEF DESCRIPTION OF DRAWINGS

FIG. 1 is a conceptual illustration of a communication system in accordance with an embodiment of the invention.

FIG. 2 is a conceptual illustration of a transmitter in accordance with an embodiment of the invention.

FIG. 3 is a conceptual illustration of a receiver in accordance with an embodiment of the invention.

FIG. 4a is a conceptual illustration of the joint capacity of a channel.

FIG. 4b is a conceptual illustration of the parallel decoding capacity of a channel.

FIG. 5 is a flow chart showing a process for obtaining a constellation optimized for capacity for use in a communication system having a fixed code rate and modulation scheme in accordance with an embodiment of the invention.

FIG. 6a is a chart showing a comparison of Gaussian capacity and PD capacity over the AWGN channel for traditional PAM-2,4,8,16,32.

FIG. 6b is a chart showing a comparison between Gaussian capacity and joint capacity over the AWGN channel for traditional PAM-2,4,8,16,32.

FIG. 7 is a chart showing the SNR gap to Gaussian capacity for the PD capacity and joint capacity over the AWGN channel of traditional PAM-2,4,8,16,32 constellations.

FIG. 8a is a chart comparing the SNR gap to Gaussian capacity of the PD capacity for traditional and optimized PAM-2,4,8,16,32 constellations, over the AWGN channel.

FIG. 8b is a chart comparing the SNR gap to Gaussian capacity of the joint capacity for traditional and optimized PAM-2,4,8,16,32 constellations, over the AWGN channel.

FIG. 9 is a chart showing Frame Error Rate performance of traditional and PD capacity optimized PAM-32 constellations in simulations involving several different length LDPC codes, over the AWGN channel.

FIGS. 10a-10d are locus plots showing the location of constellation points of a PAM-4 constellation optimized for PD capacity and joint capacity over the AWGN channel versus user bit rate per dimension and versus SNR.

FIGS. 11a and 11b are design tables of PD capacity and joint capacity optimized PAM-4 constellations for the AWGN channel in accordance with embodiments of the invention.

FIGS. 12a-12d are locus plots showing the location of constellation points of a PAM-8 constellation optimized for PD capacity and joint capacity over the AWGN channel versus user bit rate per dimension and versus SNR.

FIGS. 13a and 13b are design tables of PD capacity and joint capacity optimized PAM-8 constellations for the AWGN channel in accordance with embodiments of the invention.

FIGS. 14a-14d are locus plots showing the location of constellation points of a PAM-16 constellation optimized for PD capacity and joint capacity over the AWGN channel versus user bit rate per dimension and versus SNR.

FIGS. 15a and 15b are design tables of PD capacity and joint capacity optimized PAM-16 constellations for the AWGN channel in accordance with embodiments of the invention.

FIGS. 16a-16d are locus plots showing the location of constellation points of a PAM-32 constellation optimized for PD capacity and joint capacity over the AWGN channel versus user bit rate per dimension and versus SNR.

FIGS. 17a and 17b are design tables of PD capacity and joint capacity optimized PAM-32 constellations for the AWGN channel in accordance with embodiments of the invention.

FIG. 18 is a chart showing the SNR gap to Gaussian capacity for traditional and capacity optimized PSK constellations over the AWGN channel.

FIG. 19 is a chart showing the location of constellation points of PD capacity optimized PSK-32 constellations over the AWGN channel.

FIG. 20 is a series of PSK-32 constellations optimized for PD capacity over the AWGN channel at different SNRs in accordance with embodiments of the invention.

FIG. 21 illustrates a QAM-64 constructed from orthogonal Cartesian product of two PD optimized PAM-8 constellations in accordance with an embodiment of the invention.

FIGS. 22a and 22b are locus plots showing the location of constellation points of a PAM-4 constellation optimized for PD capacity over a Raleigh fading channel versus user bit rate per dimension and versus SNR.

FIGS. 23a and 23b are locus plots showing the location of constellation points of a PAM-8 constellation optimized for PD capacity over a Raleigh fading channel versus user bit rate per dimension and versus SNR.

FIGS. 24a and 24b are locus plots showing the location of constellation points of a PAM-16 constellation optimized for PD capacity over a Raleigh fading channel versus user bit rate per dimension and versus SNR.

FIGS. 25-27 are tables showing the performance of geometric PAM-8 constellations optimized for Joint Capacity over a Raleigh fading channel at specific SNRs in accordance with embodiments of the invention.

FIGS. 28-30 are tables listing the constellation points corresponding to the geometric PAM-8 constellation designs optimized for Joint Capacity over a Raleigh fading channel at specific SNRs listed in FIGS. 25-27.

FIGS. 31-32 are tables showing maximum ranges for the geometric PAM-8 constellation designs optimized for Joint Capacity over a Raleigh fading channel at specific SNRs listed in FIGS. 25-27.

FIGS. 33-37 are tables showing the performance of geometric PAM-8 constellations optimized for PD Capacity over a Raleigh fading channel at specific SNRs in accordance with embodiments of the invention.

FIGS. 38-43 are tables listing the constellation points corresponding to the geometric PAM-8 constellation designs optimized for PD Capacity over a Raleigh fading channel at specific SNRs listed in FIGS. 33-37.

FIGS. 44-47 are tables showing maximum ranges for the geometric PAM-8 constellation designs optimized for PD Capacity over a Raleigh fading channel at specific SNRs listed in FIGS. 33-37.

FIGS. 48-49 are tables showing the performance of geometric PAM-16 constellations optimized for Joint Capacity over a Raleigh fading channel at specific SNRs in accordance with embodiments of the invention.

FIGS. 50-53 are tables listing the constellation points corresponding to the geometric PAM-16 constellation designs optimized for Joint Capacity over a Raleigh fading channel at specific SNRs listed in FIGS. 48-49.

FIGS. 54-55 are tables showing maximum ranges for the geometric PAM-16 constellation designs optimized for Joint Capacity over a Raleigh fading channel at specific SNRs listed in FIGS. 48-49.

FIGS. 56-59 are tables showing the performance of geometric PAM-16 constellations optimized for PD Capacity over a fa Raleigh fading channel at specific SNRs in accordance with embodiments of the invention.

FIGS. 60-67 are tables listing the constellation points corresponding to the geometric PAM-16 constellation designs optimized for PD Capacity over a Raleigh fading channel at specific SNRs listed in FIGS. 56-59.

FIGS. 68-71 are tables showing maximum ranges for the geometric PAM-16 constellation designs optimized for PD Capacity over a Raleigh fading channel at specific SNRs listed in FIGS. 56-59.

FIGS. 72-73 are tables showing the performance of geometric PAM-32 constellations optimized for Joint Capacity over a Raleigh fading channel at specific SNRs in accordance with embodiments of the invention.

FIGS. 74-81 are tables listing the constellation points corresponding to the geometric PAM-32 constellation designs optimized for Joint Capacity over a Raleigh fading channel at specific SNRs listed in FIGS. 72-73.

FIGS. 82-83 are tables showing maximum ranges for the geometric PAM-32 constellation designs optimized for Joint Capacity over a Raleigh fading channel at specific SNRs listed in FIGS. 72-73.

FIGS. 84-85 are tables showing the performance of geometric PAM-32 constellations optimized for PD Capacity over a Raleigh fading channel at specific SNRs in accordance with embodiments of the invention.

FIGS. 86-93 are tables listing the constellation points corresponding to the geometric PAM-32 constellation designs optimized for PD Capacity over a Raleigh fading channel at specific SNRs listed in FIGS. 84-85.

FIGS. 94-95 are tables showing maximum ranges for the geometric PAM-32 constellation designs optimized for PD Capacity over a Raleigh fading channel at specific SNRs listed in FIGS. 84-85.

DETAILED DESCRIPTION OF THE INVENTION

Turning now to the detailed description of the invention, communication systems in accordance with embodiments of the invention are described that use signal constellations, which have unequally spaced (i.e. ‘geometrically’ shaped) points. In many embodiments, the communication systems use specific geometric constellations that are capacity optimized at a specific SNR. In addition, ranges within which the constellation points of a capacity optimized constellation can be perturbed and are still likely to achieve a given percentage of the optimal capacity increase compared to a constellation that maximizes d_(min), are also described. Capacity measures that are used in the selection of the location of constellation points include, but are not limited to, parallel decode (PD) capacity and joint capacity.

In many embodiments, the communication systems utilize capacity approaching codes including, but not limited to, LDPC and Turbo codes. As is discussed further below, direct optimization of the constellation points of a communication system utilizing a capacity approaching channel code, can yield different constellations depending on the SNR for which they are optimized. Therefore, the same constellation is unlikely to achieve the same gains applied across all code rates; that is, the same constellation will not enable the best possible performance across all rates. In many instances, a constellation at one code rate can achieve gains that cannot be achieved at another code rate. Processes for selecting capacity optimized constellations to achieve increased gains based upon a specific coding rate in accordance with embodiments of the invention are described below. Constellations points for geometric PAM-8, PAM-16, and PAM-32 constellations that are optimized for joint capacity or PD capacity at specific SNRs are also provided. Additional geometric PAM-8, PAM-16, and PAM-32 constellations that are probabilistically likely to provide performance gains compared to constellations that maximize d_(min), which were identified by perturbing the constellation points of geometric PAM-8, PAM-16, and PAM-32 constellations optimized for joint capacity or PD capacity, are also described. The constellations are described as being probabilistically likely to provide performance gains, because all possible constellations within the ranges have not been exhaustively searched. Within each disclosed range, a large number of constellations were selected at random, and it was verified that all the selected constellations provided a gain that exceeds the given percentage of the optimal capacity increase achieved by the optimized constellations relative to a constellation that maximizes d_(min) (i.e. a PAM equally spaced constellation). In this way, ranges that are probabilistically likely to provide a performance gain that is at least a predetermined percentage of the optimal increase in capacity can be identified and a specific geometric constellation can be compared against the ranges as a guide to the increase in capacity that is likely to be achieved. In a number of embodiments, the communication systems can adapt the location of points in a constellation in response to channel conditions, changes in code rate and/or to change the target user data rate.

Communication Systems

A communication system in accordance with an embodiment of the invention is shown in FIG. 1. The communication system 10 includes a source 12 that provides user bits to a transmitter 14. The transmitter transmits symbols over a channel to a receiver 16 using a predetermined modulation scheme. The receiver uses knowledge of the modulation scheme, to decode the signal received from the transmitter. The decoded bits are provided to a sink device that is connected to the receiver.

A transmitter in accordance with an embodiment of the invention is shown in FIG. 2. The transmitter 14 includes a coder 20 that receives user bits from a source and encodes the bits in accordance with a predetermined coding scheme. In a number of embodiments, a capacity approaching code such as a turbo code or a LDPC code is used. In other embodiments, other coding schemes can be used to providing a coding gain within the communication system. A mapper 22 is connected to the coder. The mapper maps the bits output by the coder to a symbol within a geometrically distributed signal constellation stored within the mapper. The mapper provides the symbols to a modulator 24, which modulates the symbols for transmission via the channel.

A receiver in accordance with an embodiment of the invention is illustrated in FIG. 3. The receiver 16 includes a demodulator 30 that demodulates a signal received via the channel to obtain symbol or bit likelihoods. The demapper uses knowledge of the geometrically shaped symbol constellation used by the transmitter to determine these likelihoods. The demapper 32 provides the likelihoods to a decoder 34 that decodes the encoded bit stream to provide a sequence of received bits to a sink.

Geometrically Shaped Constellations

Transmitters and receivers in accordance with embodiments of the invention utilize geometrically shaped symbol constellations. In several embodiments, a geometrically shaped symbol constellation is used that optimizes the capacity of the constellation. In many embodiments, geometrically shaped symbol constellations, which include constellation points within predetermined ranges of the constellation points of a capacity optimized constellation, and that provide improved capacity compared to constellations that maximize are used. Various geometrically shaped symbol constellations that can be used in accordance with embodiments of the invention, techniques for deriving geometrically shaped symbol constellations are described below.

Selection of a Geometrically Shaped Constellation

Selection of a geometrically shaped constellation for use in a communication system in accordance with an embodiment of the invention can depend upon a variety of factors including whether the code rate is fixed. In many embodiments, a geometrically shaped constellation is used to replace a conventional constellation (i.e. a constellation maximized for d_(min)) in the mapper of transmitters and the demapper of receivers within a communication system. Upgrading a communication system involves selection of a constellation and in many instances the upgrade can be achieved via a simple firmware upgrade. In other embodiments, a geometrically shaped constellation is selected in conjunction with a code rate to meet specific performance requirements, which can for example include such factors as a specified bit rate, a maximum transmit power. Processes for selecting a geometric constellation when upgrading existing communication systems and when designing new communication systems are discussed further below.

Upgrading Existing Communication Systems

A geometrically shaped constellation that provides a capacity, which is greater than the capacity of a constellation maximized for d_(min), can be used in place of a conventional constellation in a communication system in accordance with embodiments of the invention. In many instances, the substitution of the geometrically shaped constellation can be achieved by a firmware or software upgrade of the transmitters and receivers within the communication system. Not all geometrically shaped constellations have greater capacity than that of a constellation maximized for d_(min). One approach to selecting a geometrically shaped constellation having a greater capacity than that of a constellation maximized for d_(min) is to optimize the shape of the constellation with respect to a measure of the capacity of the constellation for a given SNR. Another approach is to select a constellation from a range that is probabilistically likely to yield a constellation having at least a predetermined percentage of the optimal capacity increase. Such an approach can prove useful in circumstances, for example, where an optimized constellation is unable to be implemented. Capacity measures that can be used in the optimization process can include, but are not limited to, joint capacity or parallel decoding capacity.

Joint Capacity and Parallel Decoding Capacity

A constellation can be parameterized by the total number of constellation points, M, and the number of real dimensions, N_(dim). In systems where there are no belief propagation iterations between the decoder and the constellation bit demapper, the constellation demapper can be thought of as part of the channel. A diagram conceptually illustrating the portions of a communication system that can be considered part of the channel for the purpose of determining PD capacity is shown in FIG. 4a . The portions of the communication system that are considered part of the channel are indicated by the ghost line 40. The capacity of the channel defined as such is the parallel decoding (PD) capacity, given by:

$C_{PD} = {\sum\limits_{i = 0}^{l - 1}\; {I\left( {X_{i};Y} \right)}}$

where X_(i) is the ith bit of the l-bits transmitted symbol, and Y is the received symbol, and I(A;B) denotes the mutual information between random variables A and B.

Expressed another way, the PD capacity of a channel can be viewed in terms of the mutual information between the output bits of the encoder (such as an LDPC encoder) at the transmitter and the likelihoods computed by the demapper at the receiver. The PD capacity is influenced by both the placement of points within the constellation and by the labeling assignments.

With belief propagation iterations between the demapper and the decoder, the demapper can no longer be viewed as part of the channel, and the joint capacity of the constellation becomes the tightest known bound on the system performance. A diagram conceptually illustrating the portions of a communication system that are considered part of the channel for the purpose of determining the joint capacity of a constellation is shown in FIG. 4b . The portions of the communication system that are considered part of the channel are indicated by the ghost line 42. The joint capacity of the channel is given by:

C _(JOINT) =I(X;Y)

Joint capacity is a description of the achievable capacity between the input of the mapper on the transmit side of the link and the output of the channel (including for example AWGN and Fading channels). Practical systems must often ‘demap’ channel observations prior to decoding. In general, the step causes some loss of capacity. In fact it can be proven that C_(G)≧C_(JOINT)≧C_(PD). That is, C_(JOINT) upper bounds the capacity achievable by C_(PD). The methods of the present invention are motivated by considering the fact that practical limits to a given communication system capacity are limited by C_(JOINT) and C_(PD). In several embodiments of the invention, geometrically shaped constellations are selected that maximize these measures.

Selecting a Constellation Having an Optimal Capacity

Geometrically shaped constellations in accordance with embodiments of the invention can be designed to optimize capacity measures including, but not limited to PD capacity or joint capacity. A process for selecting the points, and potentially the labeling, of a geometrically shaped constellation for use in a communication system having a fixed code rate in accordance with an embodiment of the invention is shown in FIG. 5. The process 50 commences with the selection (52) of an appropriate constellation size M and a desired capacity per dimension In the illustrated embodiment, the process involves a check (52) to ensure that the constellation size can support the desired capacity. In the event that the constellation size could support the desired capacity, then the process iteratively optimizes the M-ary constellation for the specified capacity. Optimizing a constellation for a specified capacity often involves an iterative process, because the optimal constellation depends upon the SNR at which the communication system operates. The SNR for the optimal constellation to give a required capacity is not known a priori. Throughout the description of the present invention SNR is defined as the ratio of the average constellation energy per dimension to the average noise energy per dimension. In most cases the capacity can be set to equal the target user bit rate per symbol per dimension. In some cases adding some implementation margin on top of the target user bit rate could result in a practical system that can provide the required user rate at a lower SNR. The margin is code dependent. The following procedure could be used to determine the target capacity that includes some margin on top of the user rate. First, the code (e.g. LDPC or Turbo) can be simulated in conjunction with a conventional equally spaced constellation. Second, from the simulation results the actual SNR of operation at the required error rate can be found. Third, the capacity of the conventional constellation at that SNR can be computed. Finally, a geometrically shaped constellation can be optimized for that capacity.

In the illustrated embodiment, the iterative optimization loop involves selecting an initial estimate of the SNR at which the system is likely to operate (i.e. SNR_(in)). In several embodiments the initial estimate is the SNR required using a conventional constellation. In other embodiments, other techniques can be used for selecting the initial SNR. An M-ary constellation is then obtained by optimizing (56) the constellation to maximize a selected capacity measure at the initial SNR_(in) estimate. Various techniques for obtaining an optimized constellation for a given SNR estimate are discussed below.

The SNR at which the optimized M-ary constellation provides the desired capacity per dimension

(SNR_(out)) is determined (57). A determination (58) is made as to whether the SNR_(out) and SNR_(in) have converged. In the illustrated embodiment convergence is indicated by SNR_(out) equaling SNR_(in). In a number of embodiments, convergence can be determined based upon the difference between SNR_(out) and SNR_(in) being less than a predetermined threshold. When SNR_(out) and SNR_(in) have not converged, the process performs another iteration selecting SNR_(out) as the new SNR_(in) (55). When SNR_(out) and SNR_(in) have converged, the capacity measure of the constellation has been optimized. As is explained in more detail below, capacity optimized constellations at low SNRs are geometrically shaped constellations that can achieve significantly higher performance gains (measured as reduction in minimum required SNR) than constellations that maximize d_(min).

The process illustrated in FIG. 5 can maximize PD capacity or joint capacity of an M-ary constellation for a given SNR. Although the process illustrated in FIG. 5 shows selecting an M-ary constellation optimized for capacity, a similar process could be used that terminates upon generation of an M-ary constellation where the SNR gap to Gaussian capacity at a given capacity is a predetermined margin lower than the SNR gap of a conventional constellation, for example 0.5 db. Alternatively, other processes that identify M-ary constellations having capacity greater than the capacity of a conventional constellation can be used in accordance with embodiments of the invention. For example, the effect of perturbations on the constellation points of optimized constellations can be used to identify ranges in which predetermined performance improvements are probabilistically likely to be obtained. The ranges can then be used to select geometrically shaped constellations for use in a communication system. A geometrically shaped constellation in accordance with embodiments of the invention can achieve greater capacity than the capacity of a constellation that maximizes d_(min) without having the optimal capacity for the SNR range within which the communication system operates.

We note that constellations designed to maximize joint capacity may also be particularly well suited to codes with symbols over GF(q), or with multi-stage decoding. Conversely constellations optimized for PD capacity could be better suited to the more common case of codes with symbols over GF(2)

Optimizing the Capacity of an M-Ary Constellation at a Given SNR

Processes for obtaining a capacity optimized constellation often involve determining the optimum location for the points of an M-ary constellation at a given SNR. An optimization process, such as the optimization process 56 shown in FIG. 5, typically involves unconstrained or constrained non-linear optimization. Possible objective functions to be maximized are the Joint or PD capacity functions. These functions may be targeted to channels including but not limited to Additive White Gaussian Noise (AWGN) or Rayleigh fading channels. The optimization process gives the location of each constellation point identified by its symbol labeling. In the case where the objective is joint capacity, point bit labelings are irrelevant meaning that changing the bit labelings doesn't change the joint capacity as long as the set of point locations remains unchanged.

The optimization process typically finds the constellation that gives the largest PD capacity or joint capacity at a given SNR. The optimization process itself often involves an iterative numerical process that among other things considers several constellations and selects the constellation that gives the highest capacity at a given SNR. In other embodiments, the constellation that requires the least SNR to give a required PD capacity or joint capacity can also be found. This requires running the optimization process iteratively as shown in FIG. 5.

Optimization constraints on the constellation point locations may include, but are not limited to, lower and upper bounds on point location, peak to average power of the resulting constellation, and zero mean in the resulting constellation. It can be easily shown that a globally optimal constellation will have zero mean (no DC component). Explicit inclusion of a zero mean constraint helps the optimization routine to converge more rapidly. Except for cases where exhaustive search of all combinations of point locations and labelings is possible it will not necessarily always be the case that solutions are provably globally optimal. In cases where exhaustive search is possible, the solution provided by the non-linear optimizer is in fact globally optimal.

The processes described above provide examples of the manner in which a geometrically shaped constellation having an increased capacity relative to a conventional capacity can be obtained for use in a communication system having a fixed code rate and modulation scheme. The actual gains achievable using constellations that are optimized for capacity compared to conventional constellations that maximize d_(min) are considered below.

Gains Achieved By Optimized Geometrically Spaced Constellations

The ultimate theoretical capacity achievable by any communication method is thought to be the Gaussian capacity, C_(G) which is defined as:

C _(G)=½log₂(1+SNR)

Where signal-to-noise (SNR) is the ratio of expected signal power to expected noise power. The gap that remains between the capacity of a constellation and C_(G) can be considered a measure of the quality of a given constellation design.

The gap in capacity between a conventional modulation scheme in combination with a theoretically optimal coder can be observed with reference to FIGS. 6a and 6b . FIG. 6a includes a chart 60 showing a comparison between Gaussian capacity and the PD capacity of conventional PAM-2, 4, 8, 16, and 32 constellations that maximize d_(min). Gaps 62 exist between the plot of Gaussian capacity and the PD capacity of the various PAM constellations. FIG. 6b includes a chart 64 showing a comparison between Gaussian capacity and the joint capacity of conventional PAM-2, 4, 8, 16, and 32 constellations that maximize d_(min), Gaps 66 exist between the plot of Gaussian capacity and the joint capacity of the various PAM constellations. These gaps in capacity represent the extent to which conventional PAM constellations fall short of obtaining the ultimate theoretical capacity i.e. the Gaussian capacity.

In order to gain a better view of the differences between the curves shown in FIGS. 6a and 6b at points close to the Gaussian capacity, the SNR gap to Gaussian capacity for different values of capacity for each constellation are plotted in FIG. 7. It is interesting to note from the chart 70 in FIG. 7 that (unlike the joint capacity) at the same SNR, the PD capacity does not necessarily increase with the number of constellation points. As is discussed further below, this is not the case with PAM constellations optimized for PD capacity.

FIGS. 8a and 8b summarize performance of constellations for PAM-4, 8, 16, and 32 optimized for PD capacity and joint capacity (it should be noted that BPSK is the optimal PAM-2 constellation at all code rates). The constellations are optimized for PD capacity and joint capacity for different target user bits per dimension (i.e. code rates). The optimized constellations are different depending on the target user bits per dimension, and also depending on whether they have been designed to maximize the PD capacity or the joint capacity. All the PD optimized PAM constellations are labeled using a gray labeling which is not always the binary reflective gray labeling. It should be noted that not all gray labels achieve the maximum possible PD capacity even given the freedom to place the constellation points anywhere on the real line. FIG. 8a shows the SNR gap for each constellation optimized for PD capacity. FIG. 8b shows the SNR gap to Gaussian capacity for each constellation optimized for joint capacity. Again, it should be emphasized that each ‘+’ on the plot represents a different constellation.

Referring to FIG. 8a , the coding gain achieved using a constellation optimized for PD capacity can be appreciated by comparing the SNR gap at a user bit rate per dimension of 2.5 bits for PAM-32. A user bit rate per dimension of 2.5 bits for a system transmitting 5 bits per symbol constitutes a code rate of ½. At that code rate the constellation optimized for PD capacity provides an additional coding gain of approximately 1.5 dB when compared to the conventional PAM-32 constellation.

The SNR gains that can be achieved using constellations that are optimized for PD capacity can be verified through simulation. The results of a simulation conducted using a rate 1/2 LDPC code in conjunction with a conventional PAM-32 constellation and in conjunction with a PAM-32 constellation optimized for PD capacity are illustrated in FIG. 9. A chart 90 includes plots of Frame Error Rate performance of the different constellations with respect to SNR and using different length codes (i.e. k=4,096 and k=16,384). Irrespective of the code that is used, the constellation optimized for PD capacity achieves a gain of approximately 1.3 dB, which closely approaches the gain predicted from FIG. 8 a.

Capacity Optimized Pam Constellations

Using the processes outlined above, locus plots of PAM constellations optimized for capacity can be generated that show the location of points within PAM constellations versus SNR. Locus plots of PAM-4, 8, 16, and 32 constellations optimized for PD capacity and joint capacity and corresponding design tables at various typical user bit rates per dimension are illustrated in FIGS. 10a-17b . The locus plots and design tables show PAM-4, 8, 16; and 32 constellation point locations and labelings from low to high SNR corresponding to a range of low to high spectral efficiency.

In FIG. 10a , a locus plot 100 shows the location of the points of PAM-4 constellations optimized for joint capacity plotted against achieved capacity. A similar locus plot 105 showing the location of the points of joint capacity optimized PAM-4 constellations plotted against SNR is included in FIG. 10b . In FIG. 10c , the location of points for PAM-4 optimized for PD capacity is plotted against achievable capacity and in FIG. 10d the location of points for PAM-4 for PD capacity is plotted against SNR. At low SNRs, the PD capacity optimized PAM-4 constellations have only 2 unique points, while the joint capacity optimized constellations have 3. As SNR is increased, each optimization eventually provides 4 unique points. This phenomenon is explicitly described in FIG. 11a and FIG. 11b where vertical slices of FIGS. 10 ab and 10 cd are captured in tables describing some PAM-4 constellations designs of interest. The SNR slices selected represent designs that achieve capacities={0.5, 0.75, 1.0, 1.25, 1.5} bits per symbol (bps). Given that PAM-4 can provide at most log₂(4)=2 bps, these design points represent systems with information code rates R={¼, ⅜, ½, ⅝, ¾} respectively.

FIGS. 12 ab and 12 cd present locus plots of PD capacity and joint capacity optimized PAM-8 constellation points versus achievable capacity and SNR. FIGS. 13a and 13b provide slices from these plots at SNRs corresponding to achievable capacities η={0.5, 1.0, 1.5, 2.0, 2.5} bps. Each of these slices correspond to systems with code rate R=η bps/log₂(8), resulting in R={⅙, ⅓, ½, ⅔, ⅚}. As an example of the relative performance of the constellations in these tables, consider FIG. 13b which shows a PD capacity optimized PAM-8 constellation optimized for SNR=9.00 dB, or 1.5 bps. We next examine the plot provided in FIG. 8a and see that the gap of the optimized constellation to the ultimate, Gaussian, capacity (CG) is approximately 0.5 dB. At the same spectral efficiency, the gap of the traditional PAM-8 constellation is approximately 1.0 dB. The advantage of the optimized constellation is 0.5 dB for the same rate (in this case R=½). This gain can be obtained by only changing the mapper and demapper in the communication system and leaving all other blocks the same.

Similar information is presented in FIGS. 14a-14d, and 15a-15b which provide loci plots and design tables for PAM-16 PD capacity and joint capacity optimized constellations. Likewise FIGS. 16a-16d, 17a and 17b provide loci plots and design tables for PAM-32 PD capacity and joint capacity optimized constellations.

Capacity Optimized PSK Constellations Over The AWGN Channel

Traditional phase shift keyed (PSK) constellations are already quite optimal. This can be seen in the chart 180 comparing the SNR gaps of tradition PSK with capacity optimized PSK constellations shown in FIG. 18 where the gap between PD capacity and Gaussian capacity is plotted for traditional PSK-4, 8, 16, and 32 and for PD capacity optimized PSK-4, 8, 16, and 32.

The locus plot of PD optimized PSK-32 points across SNR is shown in FIG. 19, which actually characterizes all PSKs with spectral efficiency η≦5. This can be seen in FIG. 20. Note that at low SNR (0.4 dB) the optimal PSK-32 design is the same as traditional PSK-4, at SNR=8.4 dB optimal PSK-32 is the same as traditional PSK-8, at SNR=14.8 dB optimal PSK-32 is the same as traditional PSK-16, and finally at SNRs greater than 20.4 dB optimized PSK-32 is the same as traditional PSK-32. There are SNRs between these discrete points (for instance SNR=2 and 15 dB) for which optimized PSK-32 provides superior PD capacity when compared to traditional PSK constellations.

We note now that the locus of points for PD optimized PSK-32 in FIG. 19 in conjunction with the gap to Gaussian capacity curve for optimized PSK-32 in FIG. 18 implies a potential design methodology. Specifically, the designer could achieve performance equivalent or better than that enabled by traditional PSK-4,8, and 16 by using only the optimized PSK-32 in conjunction with a single tuning parameter that controlled where the constellation points should be selected from on the locus of FIG. 19. Such an approach would couple a highly rate adaptive channel code that could vary its rate, for instance, rate 4/5 to achieve and overall (code plus optimized PSK-32 modulation) spectral efficiency of 4 bits per symbol, down to 1/5 to achieve an overall spectral efficiency of 1 bit per symbol. Such an adaptive modulation and coding system could essentially perform on the optimal continuum represented by the rightmost contour of FIG. 18.

Adaptive Rate Design

In the previous example spectrally adaptive use of PSK-32 was described. Techniques similar to this can be applied for other capacity optimized constellations across the link between a transmitter and receiver. For instance, in the case where a system implements quality of service it is possible to instruct a transmitter to increase or decrease spectral efficiency on demand. In the context of the current invention a capacity optimized constellation designed precisely for the target spectral efficiency can be loaded into the transmit mapper in conjunction with a code rate selection that meets the end user rate goal. When such a modulation/code rate change occurred a message could propagated to the receiver so that the receiver, in anticipation of the change, could select a demapper/decoder configuration in order to match the new transmit-side configuration.

Conversely, the receiver could implement a quality of performance based optimized constellation/code rate pair control mechanism. Such an approach would include some form of receiver quality measure. This could be the receiver's estimate of SNR or bit error rate. Take for example the case where bit error rate was above some acceptable threshold. In this case, via a backchannel, the receiver could request that the transmitter lower the spectral efficiency of the link by swapping to an alternate capacity optimized constellation/code rate pair in the coder and mapper modules and then signaling the receiver to swap in the complementary pairing in the demapper/decoder modules.

Geometrically Shaped QAM Constellations

Quadrature amplitude modulation (QAM) constellations can be constructed by orthogonalizing PAM constellations into QAM in phase and quadrature components. Constellations constructed in this way can be attractive in many applications because they have low-complexity demappers.

In FIG. 21 we provide an example of a Quadrature Amplitude Modulation constellation constructed from a Pulse Amplitude Modulation constellation. The illustrated embodiment was constructed using a PAM-8 constellation optimized for PD capacity for the AWGN channel at user bit rate per dimension of 1.5 bits (corresponds to an SNR of 9.0 dB) (see FIG. 13b ). The label-point pairs in this PAM-8 constellation are {(000, −1.72), (001, −0.81), (010, 1.72), (011,−0.62), (100, 0.62), (101, 0.02), (110, 0.81), (111, −0.02)}. Examination of FIG. 21 shows that the QAM constellation construction is achieved by replicating a complete set of PAM-8 points in the quadrature dimension for each of the 8 PAM-8 points in the in-phase dimension. Labeling is achieved by assigning the PAM-8 labels to the LSB range on the in-phase dimension and to the MSB range on the quadrature dimension. The resulting 8×8 outer product forms a highly structured QAM-64 for which very low-complexity de-mappers can be constructed. Due to the orthogonality of the in-phase and quadrature components the capacity characteristics of the resulting QAM-64 constellation are identical to that of the PAM-8 constellation on a per-dimension basis. The same process can be applied for constellations optimized for the Raleigh fading channel.

N-Dimensional Constellation Optimization

Rather than designing constellations in 1-D (PAM for instance) and then extending to 2-D (QAM), it is possible to take direct advantage in the optimization step of the additional degree of freedom presented by an extra spatial dimension. In general it is possible to design N-dimensional constellations and associated labelings. The complexity of the optimization step grows exponentially in the number of dimensions as does the complexity of the resulting receiver de-mapper. Such constructions constitute embodiments of the invention and simply require more ‘run-time’ to produce.

Capacity Optimized Constellations for Fading Channels

Similar processes to those outlined above can be used to design capacity optimized constellations for fading channels in accordance with embodiments of the invention. The processes are essentially the same with the exception that the manner in which capacity is calculated is modified to account for the fading channel. A fading channel can be described using the following equation:

Y=a(t).X+N

where X is the transmitted signal, N is an additive white Gaussian noise signal and a(t) is the fading distribution, which is a function of time.

In, the case of a fading channel, the instantaneous SNR at the receiver changes according to a fading distribution. The fading distribution is Rayleigh and has the property that the average SNR of the system remains the same as in the case of the AWGN channel, E[X²]/E[N²]. Therefore, the capacity of the fading channel can be computed by taking the expectation of AWGN capacity, at a given average SNR, over a fading distribution of a, such as the Raleigh fading distribution, that drives the distribution of the instantaneous SNR.

Many fading channels follow a Rayleigh distribution. FIGS. 22a-24b are locus plots of PAM-4, 8, and 16 constellations that have been optimized for PD capacity on a Rayleigh fading channel. Locus plots versus user bit rate per dimension and versus SNR are provided. Similar processes can be used to obtain capacity optimized constellations that are optimized using other capacity measures, such as joint capacity, and/or using different modulation schemes.

Geometric PAM-8, PAM-16, and PAM-32 Constellations

As described above, geometric constellations can be obtained that are optimized for joint or PD capacity at specific SNRs. In addition, ranges can be specified for the constellation points of a geometric constellation that are probabilistically likely to result in geometric constellations that provide at least a predetermined performance improvement relative to a constellation that maximizes d_(min). Turning now to FIGS. 25-95, geometric PAM-8, PAM-16, and PAM-32 constellations optimized for joint and PD capacity over the Raleigh fading channel at specific SNRs are listed. The performances of the optimal constellations are compared to the performances of traditional constellations that maximize d_(min). Ranges for the constellation points are also defined at specific SNRs, where constellations having constellation points selected from within the ranges are probabilistically likely (with probability close to one) to result in at least a predetermined performance improvement at the specified SNR relative to a traditional constellation that maximizes d_(min).

The geometric constellations disclosed in FIGS. 25-95 are defined by points y(i) such that y(i)=k(x(i)+r(i))+c. Values for x(i) and bounds on r(i) are provided in FIGS. 25-95 for PAM-8, PAM-16, and PAM-32 optimized for joint and PD capacity under Rayleigh fading channel conditions. For PAM-8 0≦i≦7, PAM-16 0≦i≦15, and for PAM-32 0≦i≦31. To achieve optimal power efficiency, c should be set to zero. The constant k can be viewed as a scaling factor that only changes the power of the constellation. For example, the constellations disclosed in FIGS. 25-95 were scaled to provide constellations with an arbitrary power of 10,912. The same constellations can be represented using any other arbitrary scaling including (but not limited) to scaling the disclosed constellations to have an arbitrary power of 16. To scale the provided constellations with an arbitrary power of 10,912 to have an arbitrary power of 16, a scaling factor of 1/26.1151297144012 can be used. By way of example, constellation design 1 from the tables in FIGS. 86-93 can be scaled to a power of 16 as follows:

Constellation Design 1 Power = 10,912 Power = 16 −37.9896 −1.4547 −30.4644 −1.1665 −25.8467 −0.9897 −26.1314 −1.0006 −17.8829 −0.6848 −18.1165 −0.6937 −18.8919 −0.7234 −18.8919 −0.7234 −4.8867 −0.1871 −4.9004 −0.1876 −4.9384 −0.1891 −4.9384 −0.1891 −8.2966 −0.3177 −8.2966 −0.3177 −8.2966 −0.3177 −8.2966 −0.3177 37.9896 1.4547 30.4644 1.1665 25.8467 0.9897 26.1314 1.0006 17.884 0.6848 18.1156 0.6937 18.8918 0.7234 18.8918 0.7234 4.8863 0.1871 4.9008 0.1877 4.9384 0.1891 4.9384 0.1891 8.2966 0.3177 8.2966 0.3177 8.2966 0.3177 8.2966 0.3177

Similarly, modifying k also scales the resulting maximum ranges for constellation design 1 as shown in the tables from FIGS. 94-95 as follows:

Constellation Design 1 Percentage of Optimal Capacity Gain 5% 20% 50% 90% 100% Maximum ranges ±2.1890 ±1.9700 ±1.5960 ±0.6870 ±0 where k = 1 Maximum ranges ±8.3821 × 10⁻² ±7.5435 × 10⁻² ±6.1114 × 10⁻² ±2.6307 × 10⁻² ±0 where k = 1/26.1151297144012

By way of further example, constellation design 8 from the tables in FIGS. 86-93 can be scaled to a power of 16 as follows:

Constellation Design 8 Power = 10,912 Power = 16 −38.8088 −1.4861 −31.7370 −1.2153 −25.0630 −0.9597 −26.1892 −1.0028 −16.6675 −0.6382 −16.6894 −0.6391 −18.4668 −0.7071 −18.4668 −0.7071 −3.7549 −0.1438 −3.7549 −0.1438 −3.7549 −0.1438 −3.7549 −0.1438 −9.2115 −0.3527 −9.2115 −0.3527 −9.0480 −0.3465 −9.0888 −0.3480 38.8190 1.4865 31.7383 1.2153 25.0617 0.9597 26.1893 1.0028 16.6677 0.6382 16.6875 0.6390 18.4662 0.7071 18.4662 0.7071 3.7539 0.1437 3.7539 0.1437 3.7539 0.1437 3.7539 0.1437 9.2107 0.3527 9.2107 0.3527 9.0476 0.3465 9.0509 0.3466

Similarly, modifying k also scales the resulting maximum ranges for constellation design 8 as shown in the tables from FIGS. 94-95 as follows:

Constellation Design 8 Percentage of Optimal Capacity Gain 5% 20% 50% 90% 100% Maximum ranges ±1.9640 ±1.7680 ±1.2890 ±0.6020 ±0 where k = 1 Maximum ranges ±7.5205 × 10⁻² ±6.7700 × 10⁻² ±4.9358 × 10⁻² ±2.3052 × 10⁻² ±0 where k = 1/26.1151297144012

By way of further example, constellation design 13 from the tables in FIGS. 86-93 can be scaled to a power of 16 as follows:

Constellation Design 13 Power = 10,912 Power = 16 −39.2317 −1.5023 −32.0637 −1.2278 −24.6366 −0.9434 −26.214 −1.0038 −15.9554 −0.6110 −15.9554 −0.6110 −18.7399 −0.7176 −18.5633 −0.7108 −3.3634 −0.1288 −3.3645 −0.1288 −3.3963 −0.1301 −3.3963 −0.1301 −9.556 −0.3659 −9.556 −0.3659 −8.9797 −0.3439 −8.9852 −0.3441 39.2363 1.5024 32.0635 1.2278 24.6352 0.9433 26.2126 1.0037 15.9543 0.6109 15.9543 0.6109 18.7381 0.7175 18.5618 0.7108 3.3646 0.1288 3.3653 0.1289 3.3966 0.1301 3.3966 0.1301 9.5562 0.3659 9.5562 0.3659 8.9803 0.3439 8.9851 0.3441

Similarly, modifying k also scales the resulting maximum ranges for constellation design 13 as follows:

Constellation Design 13 Percentage of Optimal Capacity Gain 5% 20% 50% 90% 100% Maximum ranges ±1.7080 ±1.5370 ±1.2450 ±0.5190 ±0 where k = 1 Maximum ranges ±6.5403 × ±5.8855 × ±4.7674 × ±1.9874 × ±0 where k = 10⁻² 10⁻² 10⁻² 10⁻² 1/26.1151297144012

By way of further example, constellation design 19 from the tables in FIGS. 86-93 can be scaled to a power of 16 as follows:

Constellation Design 19 Power = 10,912 Power = 16 −39.1811 −1.5003 −32.19 −1.2326 −24.2565 −0.9288 −26.6863 −1.0219 −15.3045 −0.5860 −15.3045 −0.5860 −19.4031 −0.7430 −18.8423 −0.7215 −2.6518 −0.1015 −2.652 −0.1016 −3.5662 −0.1366 −3.5662 −0.1366 −10.1865 −0.3901 −10.1865 −0.3901 −8.2708 −0.3167 −8.2773 −0.3170 39.1856 1.5005 32.1917 1.2327 24.2575 0.9289 26.686 1.0219 15.3039 0.5860 15.3039 0.5860 19.4019 0.7429 18.8439 0.7216 2.651 0.1015 2.6514 0.1015 3.5662 0.1366 3.5662 0.1366 10.186 0.3900 10.186 0.3900 8.2697 0.3167 8.2746 0.3169

Similarly, modifying k also scales the resulting maximum ranges for constellation design 19 as shown in the tables from FIGS. 94-95 as follows:

Constellation Design 19 Percentage of Optimal Capacity Gain 5% 20% 50% 90% 100% Maximum ranges ±1.3880 ±1.2490 ±1.0120 ±0.4430 ±0 where k = 1 Maximum ranges ±5.3145 × ±4.7827 × ±3.8751 × ±1.6963 × ±0 where k = 10⁻² 10⁻² 10⁻² 10⁻² 1/26.1151297144012

By way of further example, constellation design 24 from the tables in FIGS. 86-93 can be scaled to a power of 16 as follows:

Constellation Design 24 Power = 10,912 Power = 16 −38.5955 −1.4779 −32.0616 −1.2277 −24.0453 −0.9207 −27.0615 −1.0362 −15.1509 −0.5802 −15.3494 −0.5878 −19.9129 −0.7625 −18.8201 −0.7207 −1.9769 −0.0757 −1.9778 −0.0757 −4.2608 −0.1632 −4.2608 −0.1632 −10.928 −0.4185 −10.9088 −0.4177 −8.0512 −0.3083 −8.0512 −0.3083 38.5954 1.4779 32.0594 1.2276 24.0456 0.9208 27.0586 1.0361 15.1515 0.5802 15.3492 0.5878 19.9127 0.7625 18.8213 0.7207 1.9789 0.0758 1.9801 0.0758 4.2617 0.1632 4.2617 0.1632 10.9264 0.4184 10.9084 0.4177 8.0509 0.3083 8.0509 0.3083

Similarly, modifying k also scales the resulting maximum ranges for constellation design 24 as shown in the tables from FIGS. 94-95 as follows:

Constellation Design 24 Percentage of Optimal Capacity Gain 5% 20% 50% 90% 100% Maximum ranges ±1.0640 ±0.9580 ±0.7760 ±0.3410 ±0 where k = 1 Maximum ranges ±4.0743 × ±3.6684 × ±2.9715 × ±1.3058 × ±0 where k = 10⁻² 10⁻² 10⁻² 10⁻² 1/26.1151297144012

By way of further example, constellation design 26 from the tables in FIGS. 86-93 can be scaled to a power of 16 as follows:

Constellation Design 26 Power = 10,912 Power = 16 −38.3887 −1.4700 −32.0467 −1.2271 −24.023 −0.9199 −27.2059 −1.0418 −15.0232 −0.5753 −15.4053 −0.5899 −20.1597 −0.7720 −18.7135 −0.7166 −1.854 −0.0710 −1.8543 −0.0710 −4.4028 −0.1686 −4.4028 −0.1686 −11.086 −0.4245 −11.0231 −0.4221 −8.0244 −0.3073 −8.0271 −0.3074 38.3284 1.4677 32.0412 1.2269 24.0205 0.9198 27.2021 1.0416 15.026 0.5754 15.4055 0.5899 20.1578 0.7719 18.7141 0.7166 1.8561 0.0711 1.8565 0.0711 4.4049 0.1687 4.4049 0.1687 11.0873 0.4246 11.0254 0.4222 8.0268 0.3074 8.0295 0.3075

Similarly, modifying k also scales the resulting maximum ranges for constellation design 26 as shown in the tables from FIGS. 94-95 as follows:

Constellation Design 26 Percentage of Optimal Capacity Gain 5% 20% 50% 90% 100% Maximum ranges ±1.0380 ±0.9340 ±0.7560 ±0.3100 ±0 where k = 1 Maximum ranges ±3.9747 × ±3.5765 × ±2.8949 × ±1.1871 × ±0 where k = 10⁻² 10⁻² 10⁻² 10⁻² 1/26.1151297144012

By way of further example, constellation design 32 from the tables in FIGS. 86-93 can be scaled to a power of 16 as follows:

Constellation Design 32 Power = 10,912 Power = 16 −37.7602 −1.4459 −32.0093 −1.2257 −24.2219 −0.9275 −27.5785 −1.0560 −14.4076 −0.5517 −15.8834 −0.6082 −20.9075 −0.8006 −18.6078 −0.7125 −1.6412 −0.0628 −1.6645 −0.0637 −4.5506 −0.1743 −4.4707 −0.1712 −11.4808 −0.4396 −10.7002 −0.4097 −7.6972 −0.2947 −7.9842 −0.3057 37.7537 1.4457 32.0068 1.2256 24.2204 0.9274 27.577 1.0560 14.4047 0.5516 15.8829 0.6082 20.906 0.8005 18.6066 0.7125 1.6435 0.0629 1.6666 0.0638 4.5532 0.1744 4.4731 0.1713 11.4808 0.4396 10.7014 0.4098 7.7004 0.2949 7.9858 0.3058

Similarly, modifying k also scales the resulting maximum ranges for constellation design 32 as shown in the tables from FIGS. 94-95 as follows:

Constellation Design 32 Percentage of Optimal Capacity Gain 5% 20% 50% 90% 100% Maximum ranges ±0.7300 ±0.6570 ±0.5320 ±0.2540 ±0 where k = 1 Maximum ranges ±2.7953 × ±2.5158 × ±2.0371 × ±0.9726 × ±0 where k = 10⁻² 10⁻² 10⁻² 10⁻² 1/26.1151297144012

By way of further example, constellation design 37 from the tables in FIGS. 86-93 can be scaled to a power of 16 as follows:

Constellation Design 37 Power = 10,912 Power = 16 −36.8223 −1.4100 −31.6346 −1.2114 −24.3912 −0.9340 −27.5952 −1.0567 −14.6073 −0.5593 −16.4909 −0.6315 −21.3627 −0.8180 −18.963 −0.7261 −1.1401 −0.0437 −2.0669 −0.0791 −5.1427 −0.1969 −4.118 −0.1577 −12.2395 −0.4687 −10.704 −0.4099 −7.3701 −0.2822 −8.5986 −0.3293 36.8223 1.4100 31.6352 1.2114 24.3917 0.9340 27.5958 1.0567 14.6073 0.5593 16.4909 0.6315 21.3631 0.8180 18.9632 0.7261 1.1397 0.0436 2.0664 0.0791 5.1424 0.1969 4.1177 0.1577 12.2394 0.4687 10.7038 0.4099 7.3697 0.2822 8.5984 0.3292

Similarly, modifying k also scales the resulting maximum ranges for constellation design 37 as shown in the tables from FIGS. 94-95 as follows:

Constellation Design 37 Percentage of Optimal Capacity Gain 5% 20% 50% 90% 100% Maximum ranges ±0.6050 ±0.5440 ±0.4410 ±0.1930 ±0 where k = 1 Maximum ranges ±2.3167 × ±2.0831 × ±1.6887 × ±0.7390 × ±0 where k = 10⁻² 10⁻² 10⁻² 10⁻² 1/26.1151297144012

By way of further example, constellation design 43 from the tables in FIGS. 86-93 can be scaled to a power of 16 as follows:

Constellation Design 43 Power = 10,912 Power = 16 −35.7806 −1.3701 −31.1542 −1.1930 −24.5029 −0.9383 −27.5041 −1.0532 −15.019 −0.5751 −17.0214 −0.6518 −21.7164 −0.8316 −19.3443 −0.7407 −0.9323 −0.0357 −2.391 −0.0916 −5.6791 −0.2175 −4.1763 −0.1599 −12.8773 −0.4931 −11.1053 −0.4252 −7.5809 −0.2903 −9.1812 −0.3516 35.7805 1.3701 31.1543 1.1930 24.5029 0.9383 27.5041 1.0532 15.019 0.5751 17.0214 0.6518 21.7164 0.8316 19.3443 0.7407 0.9322 0.0357 2.391 0.0916 5.6791 0.2175 4.1762 0.1599 12.8772 0.4931 11.1053 0.4252 7.5809 0.2903 9.1813 0.3516

Similarly, modifying k also scales the resulting maximum ranges for constellation design 43 as shown in the tables from FIGS. 94-95 as follows:

Constellation Design 43 Percentage of Optimal Capacity Gain 5% 20% 50% 90% 100% Maximum ranges ±0.3890 ±0.3500 ±0.2830 ±0.1220 ±0 where k = 1 Maximum ranges ±1.4896 × ±1.3402 × ±1.0837 × ±0.4672 × ±0 where k = 10⁻² 10⁻² 10⁻² 10⁻² 1/26.1151297144012

By way of further example, constellation design 52 from the tables in FIGS. 86-93 can be scaled to a power of 16 as follows:

Constellation Design 52 Power = 10,912 Power = 16 −34.7205 −1.3295 −30.6582 −1.1740 −24.6158 −0.9426 −27.3954 −1.0490 −15.478 −0.5927 −17.5199 −0.6709 −22.0414 −0.8440 −19.7371 −0.7558 −0.8935 −0.0342 −2.5729 −0.0985 −6.0512 −0.2317 −4.3433 −0.1663 −13.4322 −0.5143 −11.5522 −0.4424 −7.8881 −0.3021 −9.6593 −0.3699 34.719 1.3295 30.6573 1.1739 24.6154 0.9426 27.3948 1.0490 15.4781 0.5927 17.52 0.6709 22.0413 0.8440 19.7371 0.7558 0.894 0.0342 2.5734 0.0985 6.0516 0.2317 4.3438 0.1663 13.4325 0.5144 11.5525 0.4424 7.8885 0.3021 9.5697 0.3664

Similarly, modifying k also scales the resulting maximum ranges for constellation design 52 as shown in the tables from FIGS. 94-95 as follows:

Constellation Design 52 Percentage of Optimal Capacity Gain 5% 20% 50% 90% 100% Maximum ranges ±0.2550 ±0.2290 ±0.1850 ±0.0800 ±0 where k = 1 Maximum ranges ±0.9764 × ±0.8769 × ±0.7084 × ±0.3063 × ±0 where k = 10⁻² 10⁻² 10⁻² 10⁻² 1/26.1151297144012

As one can readily appreciate, modifying k also scales the resulting maximum ranges for any geometric PAM-32 constellation designs optimized for PD Capacity over a Raleigh fading channel at specific SNRs.

In addition to optimized constellations, FIGS. 25-95 specify ranges for the points of a geometric constellation, where selecting the points of a constellation from within the ranges is probabilistically likely to provide a geometric constellation having at least a predetermined performance improvement relative to a constellation that maximizes d_(min). The ranges are expressed as a maximum value for the constellation range parameter, r(i), which specifies the amount by which the point x(i) in the constellation is perturbed relative to the location of the corresponding point in the optimal constellation. A communication system using a constellation formed from constellations points selected from within the ranges specified by the maximum value (i.e. −r_(max)≦r(i)≦r_(max)) is probabilistically likely to achieve a predetermined performance improvement relative to a constellation that maximizes d_(min). The predetermined performance improvements associated with the ranges specified in FIGS. 25-95 are expressed in terms of a percentage of the increase in capacity achieved by the optimized constellation relative to a constellation that maximizes Constellations formed from constellation points selected from within the ranges are probabilistically likely to achieve an increase in capacity at least as great as the indicated percentage.

With regard to the specific tables shown in FIGS. 25-95, each table is one of three different types of table. A first set of tables shows the performance of specific geometric constellations optimized for joint capacity or PD capacity. These tables include 6 columns. The first column enumerates a design number. The second column provides the SNR at which the constellation was optimized for the design defined by the entry in the first column. The third column provides the capacity achieved by the optimized constellation (Opt. Cap) at the SNR given in the second column. The fourth column provides the capacity achieved (Std. Cap) by a traditional uniformly spaced constellation i.e. a PAM constellation that maximizes (with the same number of points as the optimized constellation and where binary reflective gray labeling is assumed) at the SNR given in the second column. The fifth column shows the gain in bits per transmission provided by the optimized constellation over a constellation that maximizes d_(min). The sixth column shows the percentage gain in capacity provided by the optimized constellation over the capacity provided by the traditional uniformly spaced constellation.

A second set of tables lists the constellation points of the designs indicated in the first set of tables. These tables contain 9 columns. The first column enumerates a design number.

The remaining 8 columns describe a constellation point x(i) enumerated by label in the second row of the table. Labels are given in decimal number format. With PAM 8 as an example, a label of 011 is given as the decimal number 3.

The third set of tables specifies maximum perturbation ranges for the capacity optimized constellations indicated in the first set of tables, where the maximum ranges correspond to a high probabilistic likelihood of at least a predetermined performance improvement relative to a constellation that maximizes . These tables contain 8 columns. The first enumerates a design number (corresponding to a design from one of the aforementioned tables). The second column provides the SNR for the design defined by the entry in the first column. The remaining 5 columns describe parameter r_(max) which is the maximum amount any point in the designed constellation may be perturbed (in either the positive or negative direction) and still retain, with probability close to unity, at least the gain noted by each column header of the joint or PD capacity as a percentage of the gain provided by the corresponding optimized point design over a traditional constellation that maximizes (all at the given SNR). Each table has a last column showing that if 100% of the gain afforded by the optimized constellation is desired, then parameter r(i) must be equal to zero (no deviation from designed points described in the point specification tables).

Labelling of Constellations Using Cyclically Rotated Binary Reflective Gray Labels

In performing optimization with respect to PD capacity, a conjecture can be made that constraining the optimization process to the subsequently described class of labelings results in no or negligible loss in PD capacity (the maximum observed loss is 0.005 bits, but in many cases there is no loss at all). Use of this labeling constraint speeds the optimization process considerably. We note that joint capacity optimization is invariant to choice of labeling. Specifically, joint capacity depends only on point locations whereas PD capacity depends on point locations and respective labelings.

The class of cyclically rotated binary reflective gray labels can be used. The following example, using constellations with cardinality 8, describes the class of cyclically rotated binary reflective gray labels. Given for example the standard gray labeling scheme for

PAM-8: 000, 001, 011, 010, 110, 111, 101, 100 Application of a cyclic rotation, one step left, yields: 001, 011, 010, 110, 111, 101, 100, 000 Application of a cyclic rotation, two steps left, yields: 011, 010, 110, 111, 101, 100, 000, 001

For a constellation with cardinality 8, cyclic rotations of 0 to 7 steps can be applied. It should be noted that within this class of labelings, some labelings perform better than others. It should also be noted that different rotations may yield labelings that are equivalent (through trivial column swapping and negation operations). In general, labelings can be expressed in different but equivalent forms through trivial operations such as column swapping and negation operations. For example the binary reflective gray labels with one step rotation:

001, 011, 010, 110, 111, 101, 100, 000 Can be shown to be equivalent to: 000, 001, 011, 111, 101, 100, 110, 010

The above equivalence can be shown by the following steps of trivial operations:

1) Negate the third column. This gives 000, 010, 011, 111, 110, 100, 101, 001 2) Swap the second and third columns. This gives 000, 001, 011, 111, 101, 100, 110, 010

The two labelings are considered equivalent because they yield the same PD Capacity as long as the constellation points locations are the same.

In the constellation point specifications shown in FIGS. 25-95, a labeling can be interchanged by any equivalent labeling without affecting the performance parameters. A labeling used in the specifications may not directly appear to be a cyclically rotated binary reflective gray labeling, but it can be shown to be equivalent to one or more cyclically rotated binary reflective gray labelings.

Prior Art Geometric Constellations

Geometric constellations have been specified in the prior art in attempts to achieve performance gains relative to constellations that maximize d_(min). Examples of such constellations are disclosed in Sommer and Fettweis, “Signal Shaping by Non-Uniform QAM for AWGN Channels and Applications Using Turbo Coding” ITG Conference Source and Channel Coding, p. 81-86, 2000. The specific constellations disclosed by Sommer and Fettweis for PAM-8, PAM-16, and PAM-32 are as follows:

PAM-8: −1.6630 −0.9617 −0.5298 −0.1705 0.1705 0.5298 0.9617 1.6630 PAM-16: −1.9382 −1.3714 −1.0509 −0.8079 −0.6026 −0.4185 −0.2468 −0.0816 0.0816 0.2468 0.4185 0.6026 0.8079 1.0509 1.3714 1.9382 PAM-32: −2.1970 −1.7095 −1.4462 −1.2545 −1.0991 −0.9657 −0.8471 −0.7390 −0.6386 −0.5441 −0.4540 −0.3673 −0.2832 −0.2010 −0.1201 −0.0400 0.0400 0.1201 0.2010 0.2832 0.3673 0.4540 0.5441 0.6386 0.7390 0.8471 0.9657 1.0991 1.2545 1.4462 1.7095 2.1970

Another class of geometric constellations is disclosed in Long et al., “Approaching the AWGN Channel Capacity without Active Shaping” Proceedings of International Symposium on Information Theory, p. 374, 1997. The specific PAM-8, PAM-16, and PAM-32 constellations disclosed by Long et al. are as follows:

PAM-8: −3 −1 −1 −1 1 1 1 3 PAM-16: −4 −2 −2 −2 −2 0 0 0 0 0 0 2 2 2 2 4 PAM-32: −5 −3 −3 −3 −3 −3 −1 −1 −1 −1 −1 −1 −1 −1 −1 −1 1 1 1 1 1 1 1 1 1 1 3 3 3 3 3 5

The above prior art constellations are geometric and can provide performance improvements at some SNRs relative to constellations that maximize . The performance of the constellations varies with SNR and at certain SNRs the constellations are proximate to capacity optimized constellations. Therefore, the ranges specified in FIGS. 25-95 are defined so that prior art constellations are excluded at the specific SNRs at which these constellations are proximate to a capacity optimized constellation.

Constructing Multidimensional Constellations

The tables shown in FIGS. 25-95 can be used to identify optimal N-dimensional constellations. The optimized multi-dimensional constellation can be determined by finding the Cartesian power X^(n) and the resulting labeling constructed by finding the corresponding Cartesian power of L^(n). Ranges within which the multi-dimensional constellation points can be selected (i.e. perturbed), can then be defined with respect to each constellation point of the constructed multi-dimensional constellation, using an n-dimensional perturbation vector, such that each component of the perturbation vector has a magnitude that is less than r_(max) defined by the range tables.

Example of a QAM Constellation

The optimized constellation points for a PAM-8 constellation optimized for PD capacity at SNR=9 dB are as follows:

−7.3992 −4.8128 −1.3438 −2.0696 7.3991 4.8129 1.3438 2.0691

The labelings corresponding to the above PAM-8 constellation points are:

000 001 010 011 100 101 110 111

Using this PAM-8 constellation, it is possible to construct a QAM-64 constellation. While PAM-8 maps 3 bits to one dimension, QAM-64 maps 6 bits to two dimensions. The first three bits will determine the location in the X-dimension and the second three bits will determine the location in the Y-dimension. The resulting QAM-64 constellation for example will map the bits 000 010 to the two dimensional constellation point (−7.3992, −1.3438), and 111 110 to the two dimensional constellation point (2.0691, 1.3438). The points corresponding to the remaining labels can be derived in a similar manner.

The ranges shown in FIGS. 25-95 can be utilized to select QAM constellations in a similar manner to that outlined above with respect to the selection of a PAM-8 constellation based upon ranges specified with respect to a PAM-8 constellation optimized for PD capacity at 9 dB. A range of 0.282 can be applied to every component of each two dimensional constellation point (based on the 50% of the achievable gain in capacity column). For example, the two points (−7.4992, −1.1438) and (2.2691, 1.1438) are within the ranges as they are spaced distances (−0.1, 0.2) and (0.2, −0.2) respectively from the optimized constellation points. In this way, the ranges can be used to identify constellations that are probabilistically likely to result in a performance improvement relative to a constellation that maximizes d_(min).

The same procedure can apply to a constellation optimized for joint capacity. However, the choice of labeling does not affect joint capacity. The above procedure can similarly be applied to an N-dimensional constellation constructed from a PAM constellation.

Although the present invention has been described in certain specific embodiments, many additional modifications and variations would be apparent to those skilled in the art. It is therefore to be understood that the present invention may be practiced otherwise than specifically described, without departing from the scope and spirit of the present invention. Thus, embodiments of the present invention should be considered in all respects as illustrative and not restrictive. 

What is claimed is:
 1. A digital communication system, comprising: a transmitter configured to transmit signals to a receiver via a communication channel; wherein the transmitter, comprises: a coder configured to receive user bits and output encoded bits at an expanded output encoded bit rate; a mapper configured to map encoded bits to symbols in a symbol constellation; a modulator configured to generate a signal for transmission via the communication channel using symbols generated by the mapper; wherein the receiver, comprises: a demodulator configured to demodulate the received signal via the communication channel; a demapper configured to estimate likelihoods from the demodulated signal; a decoder that is configured to estimate decoded bits from the likelihoods generated by the demapper; and wherein the symbol constellation is a PAM-8 symbol constellation having constellation points within at least one of the ranges specified in FIGS. 25-47. 